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Abstract 

We study entanglement dynamics in quantum Brownian motion (QBM) models. Our 
main tool is the Wigner function propagator. Time evolution in the Wigner picture is 
physically intuitive and it leads to a simple derivation of a master equation for any number 
of system harmonic oscillators and spectral density of the environment. It also provides 
generalized uncertainty relations, valid for any initial state, that allow a characterization 
of the environment in terms of the modifications it causes to the system's dynamics. In 
particular, the uncertainty relations are very informative about the entanglement dynamics 
qh" of Gaussian states, and to a lesser extent for other families of states. For concreteness, we 

apply these techniques to a bipartite QBM model, describing the processes of entanglement 
creation, disentanglement, and decoherence at all temperatures and time scales. 
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v/-) ; 1 Introduction 

The study of quantum entanglement is both of practical and theoretical significance: entangle- 
ment is viewed as a physical resource for quantum-information processing and it constitutes a 
major issue in the foundations of quantum theory. The quantification of entanglement is difficult 
in multipartite systems (see, for example, Refs. [HEJElll]); however, there are useful separability 
criteria and entanglement measures for bipartite states, pure and mixed [5] l6l \7\ IH[ |9[ \T0 \ [TT] [T2] ) . 

Realistic quantum systems, including multipartite ones, cannot avoid interactions with their 
environments, which can degrade their quantum coherence and entanglement. Thus quantum 
decoherence and disentanglement are obstacles to quantum-information processing |13|. [141 \TE\ 
I16j . On the other hand, some environments act as intermediates that generate entanglement in 
multipartite systems, even if the components do not interact directly [171 I18|. [T9] . The theoretical 
study of entanglement dynamics in open quantum systems has uncovered important physical 
effects, such as the sudden death of entanglement [20l [21], entanglement revival after sudden 
death [22], the significance of non-Markovian effects |23[ l2"^[ [25], the possibility of a rich phase 
structure for the asymptotic behavior of entanglement [26[ [27] , and intricacies in the evolution 
of entanglement in multipartite systems [28J . 

Here, we study entanglement and decoherence in quantum Brownian motion (QBM) models 
, focusing on their description in terms of generalized uncertainty relations. Our main 
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tool in this study is the Wigner function propagator. QBM models are defined by a quadratic 
total Hamiltonian, and they are characterized by a Gaussian propagator. This propagator is 
solely determined by two matrices: one corresponding to the classical dissipative equations of 
motion and one containing the effect of environment-induced diffusion. In Sec. II we provide 
explicit formulas for their determination. 

The simplicity of time evolution in the Wigner picture leads to a concise derivation of an 
exact master equation for general QBM models, with any number of system oscillators and 
spectral density. Moreover, time evolution in the Wigner picture allows for a derivation of 
generalized uncertainty relations, valid for any initial state, that incorporate the influence of 
the environment upon the system. These uncertainty relations generalize the ones of Ref. |32j 
to QBM models with an arbitrary number of system oscillators — see also Refs. |33(. 125]. Their 
most important feature is that the lower bound is independent of the initial state, and for this 
reason, they allow for general statements about the process of decoherence and thermalization. 

The uncertainty relations are also related to separability criteria for bipartite systems [12] . 
Hence, they provide an important tool for the study of entanglement dynamics. For Gaussian 
states, in particular, the uncertainty relations, derived here, provide a general characterization 
of processes such as entanglement creation and disentanglement without the need to specify 
detailed properties of the initial state. However, uncertainty relations do not suffice to distinguish 
all entangled non-Gaussian states. For such states, the description of entanglement dynamics 
from the uncertainty relations is rather partial, but still leads to nontrivial results. 

The uncertainty relations derived in this article apply to any open quantum system charac- 
terized by Gaussian propagation, and they are expressed solely in terms of the coefficients of the 
Wigner function propagator. They can be used for the study of entanglement dynamics, not only 
in bipartite but also in multipartite systems. To demonstrate their usefulness, we apply them to 
a concrete bipartite QBM model system that has been studied by Paz and Roncanglia [261 EZ] • 
In this model, there exist two coupled subalgebras of observables, only one of which couples 
directly to the environment. For a special case of the system parameters, considered in Ref. 
[26] . one of the subalgebras is completely decoupled, and thus there exists a decoherence-free 
subspace for the system. Here we focus on the generic case, also explored in Ref. [27]. 

We find that in the high-temperature regime, decoherence and disentanglement are generic 
and the uncertainty relations allow for an identification of the characteristic timescales, which 
in some cases may be of very different orders of magnitude. At low temperature, entanglement 
creation often occurs and we demonstrate that it is accompanied by "entanglement oscillations" , 
that is, a sequence of entanglement sudden death and revivals at early times. In this regime, 
there is no decoherence, and disentanglement arises because of relaxation. At a time scale of 
the order of relaxation time the system tends to a unique asymptotic state, which coincides 
with a thermal state at the weak-coupling limit. The generalized uncertainty relations allow for 
the determination of upper limits to disentanglement time with respect to all Gaussian initial 
states. 

The structure of the article is the following. In Sec. II we construct the Wigner function 
propagator for the most general QBM model and we provide explicit formulas for the propaga- 
tor's coefficients. The master equation is then simply obtained from the propagator. In Sec. Ill 
we construct the generalized uncertainty relations valid for all QBM models, we show that they 
can be used for the study of multipartite entanglement, and we then consider their special case 
in the model of Refs. [261 127] , In Sec. IV we employ the uncertainty relations for the study of 
decoherence, disentanglement, and entanglement creation in different regimes and time scales of 



this model. 

2 Quantum Brownian motion models for multipartite systems 

In this section, we consider the most general setup for quantum Brownian motion, namely, a 
system of N harmonic oscillators of masses M r and frequencies r interacting with a heat bath. 
The heat bath is modeled by a set of harmonic oscillators of masses rrii and frequencies tOj, 
initially at a thermal state of temperature T. The Hamiltonian of the total system is a sum of 
three terms H = H sys + H env + Hint-, where 

*~ - ?(^ + ^) « 

*~ = Bs^ + ^J < 2 > 

i l 

Hint = 2-^2-^CirX r qi, (3) 

i a 

where X r and P r are the position and momentum operators for the system oscillators and qi 
and pi are the position and momentum operators for the environment oscillator. The interaction 
Hamiltonian Eq. ([3]) involves different couplings Ci r of each system oscillators to the bath. Thus 
it can also be used to describes systems different from the classic setup of Brownian motion, for 
example, particle detectors at different locations interacting with a quantum field [34]. 

For an initial state that is factorized in system and environment degrees of freedom the 
evolution of the reduced density matrix for the system variables is autonomous, and it can be 
expressed in terms of a master equation. For the issues we explore in this article, in particular 
entanglement dynamics, the determination of the propagator of the reduced density matrix is 
more important than the construction of the master equation, because it allows us to follow the 
time evolution of the relevant observables. The construction of the propagator is simpler in the 
Wigner picture. 

Instead of the density operator, we work with the Wigner function, defined by 

W(X, P) = ^p / ^Ce- iP - C p(X + ic, X - \(). (4) 

Its inverse is 

p(X,Y) = J d N P e ' p -( x - x ') W(i(X + X'),P). (5) 

For a factorized initial state, time evolution in QBM models is encoded in the density matrix 
propagator J(X/, Y/,t|Xo, Yo,0), defined by 

Pt(X f ,Y f ) = Jd N X d N Y J(X / ,Y / ,t|X ,Y ,0)po(Xo,Yo). (6) 

The Wigner function propagator is defined as 
ff(X / ,P /> t|X ,Po,0) = / ^g^e^-Co-^.0 J(X/ + |, X/ - |,t|x + |,X - |,0). (7) 



Denoting the phase-space coordinates by the vector 

£ a = (X 1 ,Pi,X 2 ,P 2 ,...,X N ,P N ), a = l,2,...,2N, (8) 

we write the Wigner function propagator compactly as -Kt(£/,£o) and express Eq. © as 

Wt(0 = J j^pmf,Zo)W (tio), (9) 

where Wt and Wq are the Wigner functions at times t and 0, respectively. 

In QBM models the Wigner function propagator is Gaussian. This follows from the fact that 
the total Hamiltonian for the system is quadratic and the initial state for the bath is Gaussian. 
The most general form of a Gaussian Wigner function propagator is 



Kt(£f,€o) = aj exp 



■k/-&(*)]<s£(t)K/ -&(*)] 



(10) 



where S~ b (t) is a positive real- valued matrix, and £d(t) is the solution of the corresponding 
classical equations of motion (including dissipation) with initial condition £ = £o at t = to- The 
equations of motion are linear, so £ c j(i) is of the form 

&(t) = Rtmi (ii) 

in terms of a matrix R^(t). 

Equation fllQj) holds if there are no "decoherence-free" subalgebras, that is, if there is no 
subalgebra of the canonical variables that remains decoupled from the environment. These 
observables evolve with a delta-function propagator, rather than with a Gaussian. However, 
this case corresponds to a set of measure zero in the space of parameters, and it can be obtained 
as a weak limit of the generic expression, Eq. (jlOp . 



In order to specify the Wigner function propagator, we must construct the matrix-valued 
functions R(t) and S(t). To this end, we consider the two-point correlation matrix V of a 
quantum state p, defined by 



Vab = \TT 



p(Ub + tea) -Tr(pi a )Tr(pi b ). (12) 



Gaussian propagation decouples the evolution of two-point correlations from any higher- 
order correlations. From Eqs. ([9]) and (I10p . we find the two-point correlation matrix, Eq. (I12h . 
Vt at time t, 

V t = R(t)V R T (t) + S(t) 1 (13) 

where Vo is the correlation matrix of the initial state. The first term in the right-hand side of 
Eq. (fT3|) corresponds to the evolution of the initial phase-space correlations according to the 
classical equations of motion. The second term incorporates the effect of environment-induced 
fluctuations and it does not depend on the initial state. Hence, the matrix S can be explicitly 
constructed, by identifying the part of the correlation matrix that does not depend on the initial 
state. 



To this end, we proceed as follows. From the Heisenberg-picture evolution of the bath 
oscillators, we obtain the equations 



&(<)+W?&(t) = X;^*r(t), (14) 



with solution 



Ut) = <fi(t) + V -^- I ds sin (cji(t - s)) X r (s), (15) 

V rriMi Jo 



where 



q?(t) = qi cos (ui{t) -\ — sin (wj£) . (16) 

rriiUJi 

For the system variables, we obtain 

2 ^ /•* - r - 



i r (t) + n2je r (t) + _ ]T / d S 7^(i-^r'(^) = E^ (i), (17) 

M r ^f JO t M r 



where 



Cj jT Ci r i 



Irr' (s) = ~ E r,^ .^ sin (^iS) ( 18 ) 



2m,; cj? 



i^j 



is the dissipation kernel. In general, the matrix 7 rr / is symmetric and has iiV(iV+l) independent 
terms, each defining a different relaxation time-scale for the system. However, symmetries of 
the couplings Cj r may reduce the number of independent components of the dissipation kernel. 
The solution of Eq. (fT7|) is 



1 1 /"* 

-MO = V(w rr /(t)X r / + — — v rT >P r i) + V — — / dsv rr '(t - s) V Cj r /g°(s), (19) 

77 W ^ M r i Jo ^ 

where v rr >(t) is the solution of the homogeneous part of Eq. (117j) . with initial conditions v rr >(0) = 
5 rr i and v rr i(0) = 0. It can be expressed as an inverse Laplace transform 

v{t) = C-\A-\z)l (20) 

where A rr /(z) = (z 2 + Q 2 )S rr i + 7g^7rr'( z ) an d lrr'{ z ) is the Laplace transform of the dissipation 
kernel. 

The classical equations of motion follow from the expectation values of X r and P r = M r X r 
in Eq. dH]) 

( X(t)\ ( v(t) v(t)M- 1 \ ( X(0) \ 

{ P(t) ) I Mv(t) Mv(t)M^ ) I P(0) } ' l ' 



where M = diag(Mi, . . . ,M r ) is the mass matrix for the system. The matrix R of Eq. (fTTI 
follows from Eq. (|2ip by a relabeling coordinated according to the definition of the vector £ a 
Eq. ©. 



We next employ Eq. (|19p . in order to construct the correlation matrix Eq. (|12p . Using 
the following equation for the correlation functions of harmonic oscillators in a thermal state at 
temperature T, 

($(')$ W))T = Sij^ coth (j±\ cos («,(« - s')) , (22) 



we find 



1 f f 

Sx r x r , = Yl U U / ds / ds'v rq (s)v qq >(s-s')v q > r >(s'), (23) 

Sp r p r , = M r M r ,Y^ / ds/ ds'v rq (s)u qq i(s-s')v q > r >(s'), ( 24 ) 

1 ft 

Sx r p r , = M r >Y^ / dsV g (s)^g'(s-s / )Vr'(s / )> (25) 

where the symmetric matrix 

v rr , (s) = ^ ^ C ^ 2 coth \-±J cos (wfa) (26) 



is the noise kernel. Similarly to the dissipation kernel, the noise kernel has ^N(N + 1) indepen- 
dent components. 

Equations (T23T [25]) together with the classical equations of motion (J2TJ) fully specify the 
Wigner function propagator. The master equation in the Wigner representation easily follows, 
by taking the time derivative of Eq. (J9]) and using the identities 

J^^-Li) a K t (C f ^o)W (Co) = -S ahd -^±i (27) 

J^ti-Zci) a (t;-tci) b m f ,to)W ({; ) = S ab + S ac S bd ^0. (28) 

The result is 

dw t_ (f?T? -i ]a d(ew t ) i. ab (f?I? -u (aocb)) d*w t (z) 

-df--( RR >o^^ + (2 s ~ {RR )cS >~deW (29) 

The method leading to the master equation (I29D is a generalization of the approach in Ref. 
|35j for the derivation of the Hu, Paz and Zhang master equation for N = 1. To the best of our 
knowledge the only other derivation of the QBM master equation in such a general setup (also 
including external force terms) is the one by Fleming, Roura and Hu, Ref. |31j. The benefit 
of the present derivation is that, by construction, it also provides the solution of the master 
equation, i.e., explicit formulas for the coefficients of the propagator. 

The first term in the right-hand side of Eq. (I29p corresponds to the Hamiltonian and dis- 



sipation terms, and the second one to diffusion with diffusion functions D (t) = (^S ~ 
RR~ 1 )cS cb '). A necessary condition for the master equation to be Markovian is that dissipa- 
tion is local, that is, that the matrix A := RR^ 1 is time independent. Then A is a generator of 
a one-parameter semi- group on the classical-state space. Moreover, the diffusion functions must 
be constant, which implies that S must be a solution of the equations S = OS + SO. 
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3 Generalized uncertainty relations 

In this section, we derive generalized uncertainty relations for the QBM models described in 
Sec. II, which are relevant to the discussion of entanglement dynamics. 

3.1 Background 

Let V, = L 2 (R N ) be the Hilbert space of a quantum system corresponding to a classical phase- 
space R . ~H carries a representation of canonical commutation relations 

[qi,Pj]=idij,i = l,...,N. (30) 

We employ a vector notation, analogous to Eq. ([SD, for the canonical operators <ji and pj. Then 
the commutation relations take the form 

[Lib]=i£lab, a,b = 1,2,... 2N, (31) 

where 



ft 



( J ... ^ 

J ... 



'--°li- <»> 



\ ... J / 

The standard uncertainty relations for this system take the form 

V > --ft. (33) 

For a bipartite system, with n degrees of freedom for the first subsystem, and N — n ones for 
the second, the Peres-Horodecki partial transpose operation defines a transformation £ — > A£, 
where A inverts the momenta of the second subsystem. Then, the correlation matrix of a 
separable state satisfies the inequality [7\ 

V > -~£l, n = AftA. (34) 



Of special interest is the case N = 2, where Eqs. (j33|) and (|34|) lead to a simple, if weaker, 
set of uncertainty relations. These have a simple generalization in the QBM model considered 
in this article. We introduce the variables 

X + = \{Xi+ X 2 ), P+ = Pi + P 2 , (35) 

X_ = ~(X!-X 2 ), P_ = P 1 -P 2 . (36) 

The partial transpose operation then interchanges P + with P_, that is, 

A(X + ,P + ,X_,P_) = (X+,P_,X_,P + ). (37) 

Hence, the uncertainty relations, 

A X+P+ := (AX + )\AP + f - F| +P+ > \, Ax_ P _ := (AX_) 2 (AP_) 2 - V^_ P _ > \, (38) 

7 



satisfied by any pair of conjugate variables (they follow from the positivity of the 2x2 diagonal 
subdeterminants of V), imply that a factorized state must satisfy the following relations 

A X+ P. := (AX+) 2 (AP_) 2 - V x+P _ > 1 Ax.p + := (AX_) 2 (AP + ) 2 - V X _ P+ > \. (39) 



If either inequality in Eq. (|39[) is violated, then the state is entangled. Hence, the uncertainty 
functions Ax + p_ and Ax_p + provide witnesses of entanglement for any state. They are weaker 
than the full Eq. (|34h . Equation (|34j) fully specifies entanglement in all Gaussian states, while 
Eq. (f39j) does so only for pure Gaussian states. 

3.2 Uncertainty relations in QBM models 

The initial correlation matrix Vq in Eq. (|13p satisfies the inequality (|33p . It follows that 



V t > -^R(t)nR T (t) + S(t). (40) 

The inequality (|40p is a generalized uncertainty relation that incorporates the effect of environment- 
induced fluctuations. It generalizes the uncertainty relations of Ref. [32J to oscillator systems 
with an arbitrary number of degrees of freedom. The right-hand side of Eq. (140p depends only 
on the coefficients of the Wigner function propagator and not on any properties of the initial 
state. Hence, Eq. (I40p provides a lower bound to the correlation matrix at time t, for a system 
that comes into contact with a heat bath at time t = 0. 

Equality in Eq. (|4T)J) is achieved for pure Gaussian states. The bound is to be understood 
in the sense of an envelope. No single Gaussian state saturates the bound in Eq. (I40p at all 
moments of time, but equality is achieved by a different family of Gaussians at each moment t. 

3.2.1 Bipartite entanglement 

When applied to a bipartite system, Eq. (fHH) implies that the condition 

- % -R{t)VLR T {t) + S{t) < -ifi (41) 

is sufficient for the existence of entangled states at time t, irrespective of the degradation caused 
by the environment. For Gaussian initial states, this condition is also necessary. 
For a factorized initial state, Eqs. (fl~3j) and (i34~|) yield 

V t > -^R(t)nR T (t) + S(t). (42) 



Inequality (j42j) is saturated for factorized pure Gaussian states, and, similarly to Eq. (J42I) . 
the lower bound to the correlation matrix is to be understood as an envelope. 

If an initially factorized state remains factorized at time t, then Vt > — 10. Then Eq. (I42p 
implies that the inequality 

- - (R{t)nR T (t) - O) + S(t) < (43) 

is a necessary condition for the preservation of factorizability at time t. 



3.2.2 Tripartite entanglement 

Equations (fT3"j) and ([30]) apply to systems of N oscillators. Used in conjunction with suitable 
separability criteria for multipartite systems [36], they also allow the derivation for uncertainty 
relations relevant multipartite systems. For example, we can use the criteria of Ref. |37| which 
apply to systems of three oscillators, labeled by the index i = 1, 2, 3. One defines the matrices 
A, that effect partial transposition with respect to the ith subsystems. Then, separable states 
satisfy 

V > ~d h tii = Afi<A, (44) 

for all i. There are some subtleties in the application of the criterion Eq. (|44p for Gaussians: 
there exist states that satisfy Eq. ([4"4"]l that are not fully separable, but only biseparable with 
respect to all possible bipartite splits — see Ref. [37] for details. However, the reasoning of Sec. 
Ill B 1 applies. 
The condition 

- i -R(t)nR T (t) + S(t)<- l -n u (45) 

for all i, is sufficient for the existence of entangled states at time t, irrespective of the degradation 
caused by the environment. For a factorized initial state, Eqs. (I13D and (I44p yield 

U t >-^(t)f^ T (t) + S(£), (46) 

for all i. Equation (|46p implies that the condition 

- l - (R{t)ttR T {t) - Oj) + S(t) < (47) 

is necessary for the preservation of factorizability at time t. 

3.3 A case model 

The uncertainty relations (|40p - (]43p hold for any Gaussian QBM system and depend only on the 
matrices R and S defining the density-matrix propagator, for which explicit expressions were 
given in Sec. II. In what follows, we elaborate on these relations in the context of a specific 
QBM model for a bipartite system, which has been studied by Paz and Roncanglia [26| I27j . 

In this model, the system consists of two harmonic oscillators with equal masses M and 
frequencies Cl\ and f^- We also consider symmetric coupling to the environment, that is., 
en = Q2 := Cj in Eq. ([3]). The latter assumption is a strong simplification, because the 
dissipation and noise kernels then become scalars, 

l(s) = Jdul(u)sm(us)( | | J, (48) 

u{s) = J duI(uj)coth(^jcos(ujs)( : j J, (49) 

where 



is the bath's spectral density. A common form for I {to) is 



. \ S 2 

to N 



I(w) =Mjto - e a*, (51) 



Co' 

where 7 is a dissipation constant, A is a high-frequency cutoff, to is a frequency scale, and the 
exponent s characterizes the infrared behavior of the bath. For this model, it is convenient to 

Q2_Q2 

employ the dimensionless parameter 5 := i J , denoting how far the system is from resonance, 
and the scaled temperature :- 



T 



In this model, the pair of oscillators is coupled to the environment only through the variables 
X + . The variable A_ is affected by the environment only through its coupling with X+, which 
is proportional to A 2 = \Q 2 — SI 2 ,]. For resonant oscillators (Oi = O2) this coupling vanishes, 
the subalgebra generated by A_ and P_ is isolated from the environment, and it is therefore 
decoherence free. This means in particular that some entanglement may persist even at late 
times. This case has been studied in detail in Ref. [26]. For nonzero A, the A_ and P_ 
subalgebra is not totally isolated from the environment. 

The uncertainty relations simplify when the environment is ohmic {s = 0). Then, dissipation 
is local and in the weak-coupling limit (7 << O,), the matrices R describing classical evolution 
take the form 

R{t) = e -^U{t), (52) 

where U{t) is a canonical transformation: U{t)£lU {t) = O. Hence, Eq. (|40p becomes 

V t >- l -e-*n + S{t). (53) 

From Eq. (j53[) we see that dissipation tends to shrink phase-space areas, but this is compensated 
by the effects of diffusion incorporated into the definition of the matrix S. For an initial factorized 
state, we obtain 

V t >~e-*F(t) + S(t), (54) 

where F{t) = U{t)ClU T {t) is an oscillating function of time. The oscillations in F{t) may lead 
to violation of the bound Vt > — §f2 for factorized states and thus to entanglement creation. 
However, the oscillating character of F{t) implies that entanglement creation will in general 
be accompanied by entanglement death and revival. For times t » 7 _1 the first term in the 
right-hand side of Eq. (f54"|) is suppressed. 

The Wigner function area. According to Eq. (|53p . the matrix Vt + §e~ 7 * — S{t) is positive. 
Its upper 2x2 submatrix in the X + , X_ coordinates should also be positive; hence, 

[(AA + ) 2 - S X+X+ }[{AP + ) 2 - S P+P+ ]} - {V X+P+ - S X+P+ ) 2 > \e- 2 ^. (55) 

By virtue of Schwartz's inequality, {AX + ) 2 S X+X+ + {AP + ) 2 S P+P+ — V x+P+ S x+P+ > 0; hence, 

Ax +P+ > \<T* + (S X+X+ S P+P+ - S 2 X+P+ ) . (56) 

10 
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Figure 1: (Color online) (a) The lower bounds for Ax + P- in Eq. (155)) and Ax + p+ in Eq. ([55)1 as 
functions of jt, for parameter values 9 = 0.7 and S = 0.38. (b) The lower bound to Ax + P- as a function 
of "ft for different values of the dimensionless temperature 0. 



Similarly, 



Ax.p. > -£* + (s x _ x _Sp_p_ 
Ax+P- > [Sx+x+Sp-P. 
•Ax- ,p+ > [Sx-X- Sp + p. 



q 2 



c-2 
D X+P- 



c2 

,5 X_P4 



(57) 
(58) 
(59) 



The uncertainty functions A Xi p correspond to the area of the projection of the Wigner 
function ellipse onto a two-dimensional subspace defined by Xi and Pj. The right-hand side of 
the inequalities are plotted in Fig. 1 as function of time. Except possibly at early times, the 
functions increase monotonically and reach a constant asymptotic value at a time scale of order 



4 Entanglement dynamics 

4.1 Disentanglement at high temperature 

A widely studied regime in quantum Brownian motion models is the so-called Fokker-Planck 
limit in ohmic environments, because in this limit the master equation is Markovian. The 
Fokker-Planck limit is defined by the condition T >> A, and then taking A — >• oo, in order to 
obtain time-local dissipation and noise. 

In this regime, thermal noise is strong, resulting in loss of quantum coherence and entan- 
glement at early times. It is convenient to work with the uncertainty functions A Xi p , because 
they can be explicitly evaluateco We find 



A x+ p + 
Ax-p- 



> i(l- 7 t + 7 2 T 2 t 4 ) 



> 



> 



[1 " 7* + 



2rp2 



rfT 



2 8 • 3 4 • 35 



A«i 12 ], 



(60) 
(61) 



J There is no loss of information in this choice, because of the rapid degradation of coherence. The sharper 
inequality, Eq. (|53[) . gives the same estimation for the characteristic time scales of these processes. 



11 



A x+ p^ > ^^AH S , (62) 



256 

7 2 T 2 



A X _ P+ > ^AV, (63) 

The above equations are obtained far from resonance for the two oscillators, that is, A >> 7. 

Equations (|60p and (|6ip represent the initial growth of fluctuations starting from purely 
quantum fluctuations at t = 0. The growth of the fluctuations for the variables X + and P+ is 
faster than that of the variables X- and P_, because the former couple indirectly to the bath. 
The — jt term in these equations indicates an initial decrease of the fluctuations, in apparent 
violation of the uncertainty principle. The violation in Eq. (|60p occurs at a timescale of order 
(7T 2 ). This is because these equations are derived taking the infinite cut-off limit A — > 00, which 
leads to violations of the positivity of the density operator at t < A -1 [32]. For t > A -1 >> T~ *, 
and T sufficiently large so that 7T 2 / A 3 >> 1, such violations do not arise. 

Ignoring the positivity-violating terms, Eq. (|60p leads to an expression tth ~ VVt^ for the 
time scale where the thermal fluctuations overcome the purely quantum ones. This is an upper 
limit to the decoherence time for the X + and P+ variables [32] . 

From Eqs. (|62p and (|63p we obtain the characteristic time scale where Ax + p_ and Ax_p + 
reach the value 2 starting from 0. This is indicative of the time scale for disentanglement t^is 
in this model: 

1 

tdis ~ ( 7 TA2)V4- (64) 

The characteristic scale for disentanglement is distinct from the time scale tth characterizing the 
growth of thermal fluctuations: 

tdis/kh = ( -Tf J • (65) 

For sufficiently small values of A, that is, weak coupling between the + and — variables, the 
disentanglement timescale may be much larger than the decoherence time scale for the X + and 
P+ variables. Hence, even if the AT_,P_ degrees of freedom are only partially protected from 
degradation from the environment, they can sustain entanglement long after the X + and P+ 
variables have decohered. 

4.2 Long-time limit 

While entanglement may be preserved much longer than the coherence of the X + and P+ degrees 
of freedom, the interaction with the environment sets the relaxation time scale 7 -1 as an upper 
limit for disentanglement time. For times t >> 7 , all states tend toward the stationary state 
Poo corresponding to a Wigner function, 



Wootf) = VdetS ™ eM-hs^Z], (66) 

where Soo is the asymptotic value of the matrix S at t — > 00. At this limit, the correlation 
matrix V coincides with S. Explicit evaluation of Eqs. (|23ti25|) shows that, as t — > 00 the only 
nonvanishing elements of the matrix S are the diagonal ones: Sx + x+,Sx-X-,Sp + p + ,Sp_p_ 
(see the Appendix). For states of this form, the uncertainty functions Ax + p_ and Ax^p + fully 
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determine entanglement. We further find that to leading order in 7/fij and fij/A, the asymptotic 
state coincides with the thermal state for Hamiltonian Hq in Eq. ([1]); hence, it is factorized. 

However, at low temperatures the thermal states are close to the boundary that separates 
factorized from entangled states (for example, they satisfy Ax + p_ — \)- Hence, the corrections 
from the nonzero values of 7/^, and Oj/A may lead the asymptotic state to retain some degree 
of entanglement, as was found in Ref. |27j . We have verified numerically that the residual 
entanglement decreases with increasing values of the cutoff parameter A. 

This result applies to a system of nondegenerate oscillations. For degenerate oscillators, the 
X_ and P_ subalgebra is protected from the environment. Hence, the asymptotic state is not 
unique and it may sustain entanglement or even be characterized by a nonterminating sequence 
of entanglement deaths and revivals. 

The analysis of Sec. II allows us to make a general characterization of the asymptotic state 
valid for any QBM model. The key observation is that the uniqueness of the asymptotic state 
is solely determined from the classical equations of motion, that is, from the matrix R^ in Eq. 
(fTTj) . In the generic case the phase space contains no dissipation- free subspace, and £^(t) — > as 
t — t- 00, irrespective of the initial condition. Hence, for times t much larger than the relaxation 
time T re i the memory of the initial state is lost from the Wigner function propagator, Eq. (jlOp . 
Moreover, if £%(t) —> sufficiently fast as t — > 00, the limit t — > 00 for the matrix S, Eqs. 
([23H25]) . is well defined. Thus a unique asymptotic state of the form ([66]) is obtained. At the 
weak-coupling limit, one expects that the asymptotic state will be close to the thermal state at 
temperature T; hence, it will be factorized. 

If, on the other hand, the classical equations of motion admit a dissipation-free subspace, 
time evolution in this subspace is Hamiltonian, and there £";(£) does not converge to a unique 
value as t — > 00. This implies that the Wigner function propagator Eq. (jlOp preserves its 
dependence on the initial variables even for t » T re \. As a consequence, an asymptotic state 
may not exist or, if it exists, it may not be unique. Hence, in this case asymptotic entanglement 
or a sequence of entanglement death and revivals is possible. 

Nonetheless, the case of a unique asymptotic state is the generic one. Dissipation-free sub- 
spaces exist only for a set of measure zero in the space of parameters (e.g., system-environment 
couplings) characterizing a QBM model. For example, even a small dependence of the coupling 
on the oscillator's position will prevent the existence of a dissipation-free subspace. Hence, 
unless some symmetry can be invoked that fully protects a subalgebra from degradation from 
the environment, we expect that the relaxation time sets an absolute upper limit to the time 
scale that entanglement can be preserved in any oscillator system interacting with a QBM-type 
environment. 

4.3 Entanglement creation 

In general, two noninteracting quantum systems may become entangled by their interaction with 
a third system. In QBM the role of the third system can be played by the environment, and 
indeed, low-temperature baths have the tendency to create entanglement. 

The uncertainty relations Eqs. (|42p and (|43p are particularly useful for the study of entan- 
glement creation. We apply them as follows. The positivity of the matrix Vt + |0 is a necessary 
criterion for a state to be factorized at time t. Hence, in a factorized state, the minimal eigen- 
value \min(t) of Vt + |fi is positive. By Eq. (j42j) . \ m in(t) is always bounded from below by the 
minimal eigenvalue of the matrix —MRQ,R T — Q) + S, which we denote as \bound(t). Hence, 
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Figure 2: (Color online) (a) The rapidly oscillating minimal eigenvalue X m i n of Vt + |fi for an initial 
factorized state |0, 1), together with the lower bound Xbound(t) corresponding to minimal eigenvalue of 
the matrix ~^(RflR T — Ci) + S, In this plot S = 0.02 and 9 = 0.21. (b) Same as in (a) but for an initial 
factorized Gaussian state. 

the function Xbmindit) determines the capacity of the environment to create entanglement irre- 
spective of the initial state. In particular, the condition that \bound{t) < implies that at least 
some factorized states can develop entanglement at time t. 

Figure 2(b) provides a plot of the minimal eigenvalue X m i n (t) of Vt + |fi for an initial 
factorized Gaussian state together with the lower bound, Xb OU nd(t), as functions of time. X m i n (t) 
oscillates rapidly at a scale of Q~ , so that at time scale of order 7 -1 we can distinguish only two 
enveloping curves that bound it from above and below. Xb oun d(t) 1S close to the lower enveloping 
curve of X m i n (t) and we note that at specific instants the inequality X m i n (t) > Xb OU nd{t) is 
saturated. 

For Gaussian states the criterion Vt < — §f2 completely specifies entanglement, hence, for 
times t that Xbound{t) > 0, no initially factorized Gaussian state can sustain entanglement. In 
Figs. 2(b) and 3, we see that Xb OU nd{t) exhibits oscillations around zero at low temperatures. This 
implies that, at least for Gaussian states, entanglement creation at low temperature is typically 
accompanied by a period of "entanglement oscillations", that is, a sequence of entanglement 
deaths and revivals, which terminates at a time scale of order 7 -1 , when the system relaxes to 
an asymptotic factorized state. 

Figure 2(a) provides a plot of X m i n (t) for an initial factorized energy eigenstate |0, 1), to- 
gether with the bound Xbound(t)- F° r non-Gaussians, a positive value of X m i n (t) does not imply 
factorizability of the state; information about entanglement is carried in higher order correlation 
functions of the system. Nonetheless, a negative value of X m i n is a definite sign of entanglement. 
Despite of the fact that X m i n (t) saturates the bound at some instants, in general its behavior is 
qualitatively different. 

In Fig. 3, the minimal eigenvalue \b OU nd(t) is plotted for different values of temperature. 
With increasing temperature the time intervals of persisting entanglement shrink and the en- 
tanglement oscillations are suppressed. At sufficiently high temperature (of order 6 > 10), no 
creation of entanglement occurs. 
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Figure 3: (Color online) The minimal eigenvalue Xbound(t) of the matrix —^(RflR —Q) + S for S = 0.02 
and different values of temperature. 

4.4 Disentanglement at low temperature 

We saw that at high temperature, the noise from the environment degrades the quantum state 
and causes rapid decoherence and disentanglement. At low temperatures (0 < 1), however, the 
noise is not sufficiently strong to cause decoherence [25], and entanglement is preserved longer. 
The physical mechanism responsible for disentanglement at low-temperature is relaxation: the 
existence of a unique asymptotic factorized state implies that at a time scale of order 7 -1 all 
memory of the initial state (including entanglement) is lost. In other words, a low temper- 
ature bath is much more efficient in creating and preserving entanglement, but relaxation to 
equilibrium will inevitably lead to a factorized state. 

By Eq. (|4*0l) . the minimal eigenvalue X m i n (t) of the matrix V t + |Q is always bounded 
from below by the minimal eigenvalue Xbound(t) 0I the matrix —^(RQR — fi) + S. Hence, the 
condition \b OU nd{t) < is sufficient for the existence of entangled states at time t. Moreover, 
the condition X m in(t) > establishes that the evolution of any Gaussian initial state at time t 
is factorized. 

Figure 4 contains plots of the minimal eigenvalue \ m in{t) of V% + |0 for two different initial 
states, together with the lower bound Xbound(t). In Fig. 4(a) the initial state is an entangled 
Gaussian, and in Fig. 4(b) the initial state is — — _. . a (|z,0) + |0, z)), where z is a coherent 

state. In both cases, \ m in(t) approaches the lower bound only after a time scale of order 7 _1 
when the system has started relaxation to a unique asymptotic state. We note that there are no 
entanglement oscillations for such states, only a gradual decay of entanglement. This behavior is 
typical for initial states that violate Eq. (I34D by a substantial margin. However, the uncertainty 
relations do not provide any significant information about the entanglement dynamics of initial 
states that are entangled, but do not violate the bound, Eq. (J34l) . This is the case, for example, 
for states of the form -^(|0, 1) + e |1, 0})). In order to study such states, we would have to 
obtain generalized uncertainty relations pertaining to correlation functions of order higher than 
2. 

In Fig. 5, we plot the minimal eigenvalue Xbound(t) as a function of time t, for different 
temperatures. As expected, the time interval during which the system sustains entangled states 
[i.e., Xbound(t) < ] shrinks with temperature. 

The uncertainty relation, Eq. (J40I) . allows for the definition of the disentanglement time tdis 
as the instant that Xb oun d{t) = 0. Thus defined, t^is is an upper bound to the disentanglement 
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Figure 4: (Color online) (a) The rapidly oscillating minimal eigenvalue of Vt + |0 for an initial 
entangled Gaussian state, together with the lower bound corresponding to minimal eigenvalue of the 
matrix -URVLR T -Q)+S. In this plot 5 = 0.02 and 9 = 0.21. (b) Same as in (a) but for the initial 
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Figure 5: (Color online) The minimal eigenvalue \bound{t) of the matrix — ^(RflR T ~ Q) + S for 5 = 0.02 
and different values of temperature. 
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Figure 6: (Color online) Disentanglement time tdis in units of 7 l as a function of the dimensionless 
temperature 9 for different values of 5. 

time for any Gaussian initial state. In general, non-Gaussian states may preserve entanglement 
for times larger than tdis- However, t<a s depends only on the matrices S and R, and the evolution 
of higher-order correlation functions of non-Gaussian states is governed by the matrices S and 
R alone. Moreover, tdis ~ 7 _1 refers to the regime of relaxation to a unique thermal equilibrium 
state, hence, the loss of any memory of the initial condition. For this reason, it is reasonable to 
assume that t^ provides a good estimation for disentanglement time that is valid for a larger 
class of initial states, at least as far as its qualitative dependence on temperature and other bath 
parameters are concerned. Figure 6 plots tdis as a function of temperature for different values 
of 5. As expected t^s decreases with temperature. However, there is no monotonic dependence 
of tdis on 5, and for 9 > 0.5, tdis is largely insensitive to 5. 

Finally, we note that the weaker uncertainty relations for the Wigner function areas Ax ± p t 
also provide an estimation for disentanglement time ti s . Since these inequalities are weaker, 
the values of tdis thus obtained are smaller, but their dependence on the parameters 5 and 9 is 
qualitatively similar. 

5 Conclusions 

The main results of our article are the following: (i) the explicit construction of the Wigner 
function propagator for QBM models with any number of system oscillators and for any spec- 
tral density; the propagator allows for a simple derivation of the corresponding master equation; 
(ii) the identification of generalized uncertainty relations valid in any QBM model that pro- 
vide a state- independent lower bound to the fluctuations induced by the environment; (hi) the 
application of the uncertainty relations to a concrete model, for the study of decoherence, dis- 
entanglement, and entanglement creation in different regimes. In particular, we showed that 
entanglement creation is often accompanied by entanglement oscillation at early times and that 
the uncertainty relations provide an upper bound to disentanglement time with respect to all 
initial Gaussian states. 

In our opinion, the most important feature of the techniques developed in this article is 
that they can be immediately generalized for addressing more complex systems and issues in 
the study of entanglement dynamics, for example, in the derivation of uncertainty relations 
for higher-order correlation functions, or for information-theoretic quantities that contain more 
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detailed information about entanglement of general initial states, and in the exploration of 
entanglement dynamics in multipartite systems and of the dependence of entanglement on the 
spatial separation of multipartite systems. 

A The coefficients in the Wigner function propagator 

In this appendix, we sketch the calculations of the coefficients in the Wigner function propagator 
for the model presented in Sec. Ill C. 

We first compute the function v(s) of Eq. (|20p in the X+, X- coordinates. To leading order 
in 7 for the poles in the Laplace transform (|2Up . we obtain 

\il 2 sin(ftis) ( 7 2 - 2il\ + 2fi| 



++VJ 4iiiii 2 (n 2 - n 2 ) 

+4 7 ^i^2 (cos(il 2 s) - cos(Ois)) - ill sm(il 2 s) (j 2 + 2il\ - 20|)1 (67) 

e -h s 

v+-(s) = [il 2 sin(Ois) - Hi sm(il 2 s)] (68) 

Z\ll\l 2 

t,_(.) = 4r ^gr_ f% ["2 Mil lS ) ( 7 2 - 87^ - 2ill + 20,3) 

+47^1^2 {cos(il 2 s) - cos(ilis)) - Hi sm(il 2 s) (^ 2 - 87^1 + 2il\ - 2fi|) (69) 

From v(s) one constructs the matrix R using Eq. (|21j) and the matrix S using Eqs. (|23H25h . 
To obtain the asymptotic state, we compute S at the limit t — > 00. The off-diagonal elements 
in the X + ,X- basis vanish, while 

Sx + x + = j^ J X ' duujf(Lo)(-2uj 2 + il\ + illf, (70) 

Sp P = ^Ml f°° duuj z f{u){-2uj 2 + il 2 + il 2 ) 2 , (71) 

7T Jo 

(■00 

(il 2 -il 2 ) 2 / duu,f(u) (72) 

JO 



Sx 



7 

M7r' 

4M 7/o2 n2 v 2 /N 



S P _ P _ = L ( n i~ n 2) / duw'fiuj), (73) 

vr Jo 



where 



^ = e-^coth(^) 

A j [2(a; 2 -^) 2 + 7 2 (u; 2 + Q2)][2( a ;2_o2)2 +7 2( u; 2 + n 2 )] - l'^ 

The asymptotic values for 5 can be evaluated to leading order in 7 /Hi and Hi/ A, by substituting 
the Lorentzians in the integrals with a delta function, that is, \{x — a) 2 + 7 2 ] -1 — j-5(x — a). 
This corresponds to the weak-damping limit of Ref . [33] . We then obtain 

1 / cothgfr cothgA 

5x + x + - 5x_x_ - ^ \—^~ + -j^-J , (75) 

Sp + p + = S p _p_ =yffii coth -^ + il 2 coth -^ J , (76) 

which correspond to an asymptotic thermal state for the pair of oscillators. 
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